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We investigate the effect of the axion cooling on the nucleosynthesis in a massive star with 16Mq 
by standard stellar evolution calculation. We find that the axion cooling suppresses the nuclear 
reactions in carbon, oxygen and silicon burning phases because of the extraction of the energy. 
As a result, larger amounts of the already synthesized neon and magnesium remain without being 
consumed to produce further heavier elements. Even in the case with the axion-photon coupling 
constant ga-/ = 10“^^ GeV“^, which is six times smaller than the current upper limit, the amount of 
neon and magnesium that remain just before the core-collapse supernova explosion is considerably 
larger than the standard value. This implies that we could give a more stringent constraint on Qa-y 
from the nucleosynthesis of heavy elements in massive stars. 


I. INTRODUCTION 

The standard model of particle physics (SM) well ex¬ 
plains general properties of the results of collider experi¬ 
ments. For example, the lattice QCD calculation, which 
is the simulation based on the first principle of the quan¬ 
tum chromodynamics (QCD), can predict the mass of 
baryons and mesons only from a few input parameters. 
On the other hand, the Lagrangian density of QCD has 
a term which violates the CP symmetry and involves the 
finite electric dipole moment (EDM) of neutron. Because 
this EDM has never been detected, it is generally con¬ 
sidered that QCD has a fine tuning problem named the 
strong CP problem. Peccei and Quinn suggested that 
the existence of an undiscovered pseudo-scalar particle 
which is associated with another C/a(1) symmetry for SM 
can solve this problem [l[ . This pseudo-scalar particle is 
named axion and have interactions with baryons, leptons 
and photons (see. e.g. [l-Q)- The coupling constant, 
ga~f, of axions to photons is related to the energy scale of 
the symmetry breaking fa as 

" 2nfa ’ 

where is a model-dependent constant. For KSVZ ii 
and DFSZ Q scenarios, \C^\ = 1.9, and 0.7, are adopted, 
respectively and several constraints on them have been 
set {e.g. [3,13])- addition, axions, which have a fi¬ 
nite mass as a results of the symmetry breaking, are a 
candidate of the cold dark matter {e.g. i). 

The interaction of axions with photons is supposed to 
affect the structure and the evolution of a star. Because 
the predicted mass of an axion is smaller than the typical 
temperature of the stellar interior, axions are expected to 
be easily produced in stars through the interaction with 
photons [2q. The conversion from photons to axions re¬ 
moves the heat in the stellar interior, which possibly gives 
an impact on the stellar structure, whereas its reaction 
rate strongly depends on temperature. 

Various possibilities concerning axions in the stellar 
interior have been explored in a wide range of stellar 
mass. By comparing the photon luminosity of the sun 


with the nuclear reaction rate that is calibrated from 
the neutrino luminosity, Gondolo and Raffelt set a con¬ 
straint ga-y < 7 X 10“^° GeV“^ Q. Tighter constraints 
can be obtained for stars in later evolutionary stages be¬ 
cause the temperature in the interior is higher than the 
temperature in main sequence stars, e.g. the sun, and 
the production rate of axions is larger as well. For ex¬ 
ample, a constraint is derived from number counts of 
horizontal branch stars; the generation of axions tends 
to shorten the duration of the horizontal branch phase, 
which contradicts to the standard stellar model without 
the effect of axions that reproduces the observed dis¬ 
tribution of horizontal branch stars within 10 % accu¬ 
racy. From this observational requirement, Ayala gives 
ga-y < 0.66 X 10“^° GeV“^0. Massive stars also give 
tight constraints on ga-y, since the interior temperature 
is suitable for the generation of axions. If one takes into 
account axions, the duration of the Helium burning is 
shorten, which would erase the blue loop stage in the 
Hertzsprung-Russell (HR) diagram required for observed 
Gepheid variable stars. By considering this effect for 
stars with 8-12 Mq, where Mq is the solar mass. Fried- 
land et al. found gay < 0.8 x 10“^° GeV“^ [l^ with MESA 
pm, a public code for one-dimensional calculations of 
stellar evolution. 

In this paper we consider the effect of axions in more 
massive stars. Pantziris & Kang estimated a constraint 
on the axion cooling rate for such massive stars by us¬ 
ing a sim ple one-zone model instead of realistic stellar 
structure [l^. In contrast, we focus on the effect of ax¬ 
ions on the nucleosynthesis of heavy nuclei at the very 
late phase of the stellar evolution just before the core¬ 
collapse supernovae. Heavy elements such as silicon, sul¬ 
fur and iron are synthesized during the last few months 
{see. P). Because of the considerably short duration of 
the nucleosynthesis of these elements, the photons gen¬ 
erated through the nuclear reactions, which do not have 
enough time to travel to the stellar surface, can be hardly 
observed. A part of the synthesized heavy nuclei are fi¬ 
nally ejected and affects elemental abundances of next 
generation stars and the chemical evolution of galaxies. 
The produced amounts of these heavy metals have the 
information of the high-temperature and high-density en- 
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FIG. 1: The Feynman diagram of the Primakoff process. Z, e 
represent ions and the electron which provide a photon via the 
magnetic field. 


vironment of the stellar interior. 

In this study, we focus on the Primakoff process [l^, 
namely axions produced by the conversion from two pho¬ 
tons (figured]), and estimate the effects of axion cooling 
on the nucleosynthesis at the late phase of massive stars 
with the MESA code. This paper is organized as follows. 
Our treatment of the axion cooling in the stellar evolu¬ 
tion code is shown in §2. In §3, we compare the structure 
of a massive star with axions to a standard case that does 
not take into account axions. In §4, we discuss the effect 
of the axion cooling on the explosive nucleosynthesis. §5 
concludes this paper. 


II. SETUP 


the region where the radiation pressure dominates, since 
the radiation pressure is proportional to T^, one can find 

M 


pi/3 \4a^l) ’ 

where TZ and a are the gas constant and the Stefan- 
Boltzmann constant, respectively. Because the value of 
the right hand side is almost stationary in the stellar 
evolution, one can regard T^/p as a constant and find 
that ea(NR&ND) oc T'^. This dependence is mentioned by 
several authors (e.g. [l^[l^.[T^b 

On the other hand, when one considers the nucleosyn¬ 
thesis of heavy elements, the temperature is so high that 
the relativistic effect needs to be taken into account, al¬ 
though the electrons are still non-degenerate [2l|, . Al- 
therr et al. reported that the formula which is valid in the 
limit of the non-relativistic and non-degenerate plasma 
can be used to the relativistic plasma, whose temperature 
exceeds the rest mass of an electron me, i.e. k^T ^ mgC^ 
[ 2 ^ In addition, as the plasma frequency wq increases, the 
emission rate suffers an exponential cut-off, oc exp(—/kuo/ 
kuT) {see. 0). Hence we adopt the following formula for 
the axion cooling rate 


£a — ea(NR&ND) ^Xp 


k^T 


( 5 ) 


= 27.2g^^,oT8V3-'i"(fcs,T)exp - 


^0 

k^T 


[erg/g/sec] 


where 30710 = 5 a 7 / 10 ~^° GeV“^. Tg and ps are tem¬ 
perature normalized by 10® K and density normalized by 
10 ® g/cm®, respectively. 


A. Axion Cooling 

The axion cooling rate per unit mass £a in hot plasma 
has been derived by several authors {e.g. [131 
They showed that the Primakoff process plays a primary 
role in the axion cooling when electrons in stars are non- 
relativistic. In particular, if the plasma frequency Wq 
is small enough to satisfy the condition, hujQ <C k^T, 
besides the non-relativistic and non-degenerate condition 
fulfilled, the emission rate can be obtained analytically 
as 

ea(NR&ND) = F{ks,T) (2) 

where ks is the Debye-Hiickel wave number which is de¬ 
fined in the Raffelt 0,[H,[23. The function F(fes,T) is 
defined as 

= ^5 i +'‘“d” (1 + ^) • 

( 3 ) 

where x = huj/k^T and k = 2ck^/hk'B,T [27j| . In horizon¬ 
tal branch stars and Oxygen burning stars, = 2.5, In 


B. Stellar Evolution 

We include the cooling by axions, eq.®, in the stellar 
evolution code MESA. We add the extra term for the ax¬ 
ion cooling in the energy transfer equation that is one of 
the basic equations governing the evolution of the stellar 
structure. The axion cooling term simply removes the 
luminosity carried by the photons emitted as a result of 
the nuclear reactions in the stellar interior. Therefore, it 
reduces the radiation pressure by these photons to mod¬ 
ify the momentum balance, and accordingly, changes the 
stellar structure if the effect is not negligible. 

We calculate the evolution of a star with M = 16Mq, 
with the solar elemental abundances from the zero-age 
main sequence phase when the hydrogen burning reac¬ 
tion is ignited at the center of the star. We take into 
account the mass loss by radiation pressure-driven stel¬ 
lar wind with an empirical mass loss rate [ 2 ^ . We follow 
the time evolution until the gravitational core-collapse 
sets in just before the supernova explosion. In addition 
to the cases with the axion cooling, we also calculate the 
standard case that does not include the effect of axions 
for comparison. 
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FIG. 2: Evolutionary tracks of stars with 16Mq that take 
into account axion coolings with different coupling constants, 
Qa-i, in an HR diagram. The long dashed-dashed (brown), 
long-dashed (red) and dashed (orange) represent the evolu¬ 
tionary track with ga-y ~ 10“^°, and 10“^^ GeV“^, 

respectively. The solid (purple) line indicates the case with¬ 
out axion cooling. 


III. RESULTS 


A. Overview of Stellar Evolution 


A star changes its luminosity and surface temperature 
with time. After the exhaustion of the hydrogen in the 
central core, a star evolves to a red giant and eventually 
the helium burning sets in to synthesize heavier elements. 
Massive stars with M > IOMq continues through the 
oxygen burning to the silicon burning phase with their 
core being non-degenerated, iron group elements domi¬ 
nate the core finally before the core-collapse supernova 
(see. H). With the stellar evolution, the star with the 
initial mass of 16 Mq lost ~ 2Mq during its lifetime by 
radiation pressure-driven stellar wind (see. (^). 

We study the effect of the axion cooling in the 
Hertzsprung-Russell (HR) diagram. In figure [21 we 
plot evolutionary tracks in cases with gaj = and 

10“^^ GeV~^ in comparison with the case without axion 
cooling [ 25 . One can hnd that these three cases show 
almost no difference. The main reason is that the du¬ 
ration of the carbon, oxygen and silicon burning phases, 
which are largely affected by the axion cooling, is sev¬ 
eral thousands years. This is much shorter than that of 
the thermal time scale of the star (~ a hundred thousand 
years.) [l^. Namely, the effect of the nuclear reaction 
deeply in the stellar interior is not still observable be¬ 
cause the travel time of the photons to the stellar surface 
is much longer. Therefore, the effect of the axion cannot 
be observed in the HR diagram. 
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FIG. 3: The abundance of alpha elements as a function of 
mass radius in units of the solar mass, Mt/Mq, at the end 
of the silicon burning phase. Long-dashed (red), dashed 
(orange), dotted (dark-green), dot-dashed (blue) lines cor¬ 
respond to the cases with Qay = 10“^^, 10“^^, 10~^®, and 
10“^^, respectively. Solid line represent the case without ax¬ 
ion cooling. 


B. Nucleosynthesis 

We investigate how the axion cooling affects the nucle¬ 
osynthesis through the evolution of the massive star. 

We compare abundances of alpha elements, oxygen 
(^®0), neon (^°Ne), magnesium (^^Mg), and silicon 
(^®Si), in addition to iron (®®Fe) of the star with dif¬ 
ferent ga-y and initial mass of M = 16 Mq in figure O 
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The abundance of ^°Ne and are enhanced for large 
ga-y near ~ “^Mq, whereas the behaviors are compli¬ 
cated. On the other hand, the abundance of ^®Si shows 
the opposite trend. The produced amount of ^®Fe is not 
influenced significantly by the axion cooling even with 
large gaj > 10“^^ GeV“^. 

In figure 01 we present the total mass of each element 
left just before the core collapse (type II) supernova, 
where each value is normalized by that of the standard 
case without axion cooling. Metals heavier than ^^Mg 
are converted to iron-group elements durin g th e explo¬ 
sive nucleosynthesis just after the explosion [25l| . Hence, 
as for the heavy elements from ^®Si to ®®Ni, we consider 
the total amount of them. 



FIG. 4: The abundance of each element i {M{Zi)) normal¬ 
ized by that from the standard case without axion cooling 
{Mo{Zi)) just before the core collapse. Long-dashed (red), 
dashed (orange), and dotted (dark-green) lines correspond to 
the cases with Qa-y = 10“^^, and GeV~^, respec¬ 

tively, where gio = sra7/(10“^°GeV“^). ’’Heavy” means the 
total amount from ^®Si to ®®Ni. 

The amount of is considerably increased for large 
gay] in the case with gay = 10“^^ the amount is enhanced 
more than three times. This trend is also weakly seen 
in the amount of ^'^Mg. On the other hand, the total 
amount of the heavier elements is smaller for larger gay 

To understand these results, we plot in figure [5] the 
time evolution of both the total energy generation rate 
and the axion emission rate at the center from the ZAMS 
to the core collapse. In the weak coupling cases with 
9ay ^ I 0 “^° GeV“^ we are considering, the axion cooling 
does not change the lifetime of the star. In such circum¬ 
stances, the effect of axion cooling is simply the extrac¬ 
tion of the thermal energy of the star, which slows down 
the nuclear reactions!^. Our calculation shows that 



FIG. 5; The time evolution of axion emission rate per unit 
mass from ZAMS (right side at tcoi — t = 10^ yr) to core¬ 
collapse (left side to tcoi — t 0) with 16 Mq, where tcoi 
is the time of the core collapse.. Long-dashed (red), dashed 
(orange), dotted (dark-green), dot-dashed (blue) lines repre¬ 
sent the cases with gay ~ 10“^^, 10“^®, and 

respectively. Solid (purple) line represents the energy genera¬ 
tion rate due to the nuclear reactions from the standard case 
without axion cooling. 

the temperature at the central region with gay = 10“^^ 
GeV“^ is 1 % lower than the standard case without ax¬ 
ion cooling at the time of the termination of the Sili¬ 
con burning at the center. In general, nuclear reaction 
rates strongly depend on temperature [e.g. 0111) and 
the small decrease of the temperature considerably slows 
down the nuclear reactions. In fact, the peak value of 
the nuclear reaction rate is suppressed by half due to the 
axion cooling at that time. 

One can find in figure [5] that the axion cooling rate 
increases monotonically over time at the central region. 
Therefore, the nuclear reactions in later phases such as 
carbon, oxygen and silicon burning phases are suppressed 
by the axion cooling and the larger amount of ^'^Ne and 
^^Mg is left, avoiding nuclear burning to synthesize heav¬ 
ier elements. Hence, the abundance of these elements, 
which are synthesized during these phases, are enhanced 
by the axion cooling as a result of the inactivation of the 
nucleosynthesis. Since these elements are not so affected 
by the later explosive nucleosynthesis [ 1 ^, the ejected 
mass would be also larger for larger gay to possibly affect 
the chemical evolution of galaxies. 

IV. CONCLUSION 

In this paper, we have studied the effect of the axion 
cooling on the nucleosynthesis in the massive star with 
I 6 M 0 by using the stellar evolution code, MESA. We have 
found that the axion cooling suppresses the nucleosynthe¬ 
sis in the carbon, oxygen and silicon burning phases even 
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in the weak coupling, 5^7 < 10“^^ GeV“^. As a result, 
the abundances of oxygen, neon and magnesium increase 
as the coupling constant ga-y increases. Even in the case 
gay = 10“^^ GeV“^, which is six times smaller than the 
current upper limit of gay = 10 “^°, the final amount of 
neon is enhanced more than three times larger than the 
standard value that does not take into account the effect 
of axion cooling. 
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Electrons, heavy ions and nucleons in stars are also ex¬ 
pected to play an important role in the axion production 
process. However we focus only on axions generated by 
photons. 

h, fcn and c are the reduced Planck constant, the Boltz¬ 
mann constant and the speed of light, respectively. 

Although Qay = 10~^° GeV“^ has already been excluded 
in a number of observations {e.g. [Til. [T2]L we plot this 
case for an illustrative purpose, in order to show the effect 
of the axion cooling. 

In general cooling increases the temperature of a star be¬ 
cause the gravo-thermal specific heat is negative. How¬ 
ever, in this case, the duration of the nuclear burning, 

0(10®) years, is much shorter than Kelvin-Helmholtz 
time scale of the star, 0(10®) years. Therefore, the ax¬ 
ion cooling of the star is too rapid to change the inner 
structure of the star, and invokes the decrease of the tem¬ 
perature of the star. 









